Dynamic properties of a diluted pyrochlore cooperative paramagnet (TbpYi_p)2Ti207 
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Investigations of the spin dynamics of the geometrically frustrated pyrochlore (TbpYi_p)2Ti207, 
using muon spin relaxation and neutron spin echo, as a function of magnetic coverage p, have 
been carried out. Our major finding is that paramagnetic fluctuations prevail as T — > for all 
values of p, and that they are sensitive to dilution, indicating a cooperative spin motion. However, 
the percolation threshold pc is not a critical point for the fluctuations. We also find that the low 
temperatures spectral density has a 1// behavior, and that dilution slows down the spin fluctuations. 



The study of geometrically frustrated magnetic sys- 
tems has introduced important new concepts to con- 
densed matterphysics including: order-by-disorder 0,0, 
spin ice [E Q ly| ; frustration driven distortion , and 
more. A recent and important addition is the concept of 
cooperative paramagnetism, whereby an interacting spin 
system dynamically fluctuates even as T ^ Q. Co- 
operative paramagnetism on the pyrochlore lattice was 
first introduced in the study of Tb2Ti2 07, where despite 
a substantial antiferromagnetic Curie- Weiss temperature 
of 20 K, a large Tb'^^ moment of 9.4 fis fluctuates contin- 
uously down to T = 50 mK. Bulk and local probes have 
been used to study the magnetic nature and the lack of a 
phase transition in Tb2Ti207 8, 9J, and the results have 
been interpreted in terms of a cooperative spin system. 
However a unique experimental identification that the 
motion of one spin in Tb2Ti207 influences the motion of 
its neighbors, is still lacking. Demonstrating conclusively 
that the fluctuations in Tb2Ti207 are indeed a result of 
spin interactions, and testing to what extent the fluctu- 
ations are cooperative, is the main aim of this work. 

Using muon spin relaxation (/iSR) and neutron spin 
echo (NSE) we investigate the dynamic fluctuations on 
a wide time range in (TbpYi„)3)2Ti207 as a function of 
p. Here the non-magnetic Y ion replaces the magnetic 
Tb ion on the corner sharing tetrahedra which make up 
the pyrochlore lattice We vary p from 1, namely, 

a pure pyrochlore, to 0.21, \yhich is below the percola- 
tion threshold Pc = 0.39 /xSR and NSE allow us 
to determine the p dependence of the waveform and pa- 
rameters for two important correlation functions (which 
are related). These are the auto correlation function 
(Si(i)Si(O)) where is the spin operator on the i'th site, 
and the intermediate scattering function S(q, t) where q 
is the wavevector. If interactions are important, we ex- 
pect the dynamic properties to be sensitive to p. More- 
over, if an infinite number of spins participate in the fluc- 
tuations we expect a dramatic change of the correlation 
functions at pc since at lower concentration the forma- 
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FIG. 1: (a) The muon asymmetry versus time at different 
temperatures for the most diluted sample. Between 102 and 
42 K and below 3 K there is hardly any temperature de- 
pendence, (b) and (c) The muon asymmetry versus time at 
T = lOOmK and various different fields for two samples. In 
(c) the error bars are omitted for clarity. The fits are to Eq.0 
as described in the text. 



tion of an infinitely large island of spins is not possible. 
Our major finding is that, indeed, the spin fluctuations 
in Tb2Ti207 result from interactions. However pc is not 
a critical point, suggesting that spin groups of finite size 
fluctuate cooperatively. We also found that the time de- 
pendent part of the correlation function is a power law. 
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Our powder samples were made by the standard solid 
state method described elsewhere Q and characterized 
by room temperature X-ray diffraction and temperature 
dependent magnetization. First we report on our fiSH 
experiments, which were conducted at the ISIS pulsed 
muon facility, Rutherford Appleton Laboratory, UK. In 
these experiments we measure the time and field de- 
pendent asymmetry A{H, t) of the decay positrons from 
muons implanted in the sample. This asymmetry is pro- 
portional to the muon polarization P{H, t) along the z 
direction, which is the direction of the external field H. 
In Fig. ^ we present the asymmetry under three differ- 
ent conditions. The log- log plot is used to emphasize the 
difference in the asymmetries between different runs. In 
panel (a) we depict the muon asymmetry versus time at 
different temperatures for the most diluted sample. Note 
that between 102 and 42 K, and below 3 K, the muon po- 
larization is almost temperature independent. In panel 
(b) and (c) we show the time dependence of the muon 
asymmetry at T = 100 mK for two different samples in 
various applied fields. At p = 0.775 the muon relaxation 
clearly decreases with increasing field. At a lower con- 
centration p = 0.21 the relaxation is lower than that at 
p = 0.775 and again decreases between 100 and lOOOG, 
but then saturates. In panel (c) the error bars are omit- 
ted for clarity. Temperature and field dependent raw 
data for the pure sample (which agree with ours) can be 
found in Ref. 8]. In all cases the raw data are fitted to 
a stretched exponential relaxation function 



A{H,t) = ^oexp - {t/Ti) 



^/3 



Br, 



(1) 



where Bg is a background signal representing muons that 
missed the sample, is the initial asymmetry, Ti is 
the spin lattice relaxation time, and (3 is the stretching 
exponent. The fits were done with a common (3 for all 
temperature and fields for a given sample. The solid 
lines in Fig. ^ S'l'e from these fits. It is important to 
mention that Ti obtained by such a fit is identical to that 
determined by the 1/e criteria {A{H,Ti) = A{H,0)/e). 

The muon relaxation rate 1 /Ti obtained from the fits 
at = 50 G as a function of temperature for several 
samples is presented in Fig. |2| The observed temper- 
ature dependence of 1/Ti is very unusual. First of all 
1/Ti shows no T dependence at the lowest T region for 
all p. This behavior is usually associated with param- 
agnetism, and it is observed in several frustrated mag- 
nets mEllIl, 

low dimensional systems with singlet 
ground state and dangling bonds [l5| . and in high spin 
molecules 16j. Most magnets show a 1/Ti maximum at 
some critical temperature. Second, the slowing down of 
spin fluctuations occurs in two steps, one at high temper- 
atures, followed by a plateau around 50 K, and another 
step starting around ~ 20 K. This might be a result of 
populations of a crystal field doublet which is 18 K above 
the ground state doublet |^. Finally, both at very high 
and very low temperatures there are significant differ- 
ences in 1/Ti between samples, however at 4 K, 1/Ti of 
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FIG. 2: Temperature dependence of the muon relaxation rate 
l/Ti at H=50G. 



all samples is the same. This is very intriguing and at 
present we have no explanation for this result. However, 
crystal field levels must be part of this explanation. In 
contrast, the p dependence of l/7i below 4 K is an in- 
dication that correlations are again playing a role in this 
low temperature region. 

In Fig. 13 we present the muon relaxation time Ti at 
T = 100 mK as a function of the longitudinal field for five 
samples. In all cases Ti increases linearly as a function of 
fields at low fields. The solid lines are linear fits to the low 
field region. This region is not affected by diluting the 
pure system by 30%, but further dilution has a strong 
impact on Ti{H). The dashed lines are continuations 
of the linear fits and serve as guides to the eye. The 
range in which Ti{H) is linear decreases with decreasing 
p, and at low concentrations Ti clearly saturates at H ^ 
0.5 kG. More importantly, the more diluted the sample, 
the greater is the sensitivity of Ti to H. This means 
that diluting the magnetic lattice slows down the spin 
fluctuation since the greatest sensitivity of Ti to external 
held is when the spin system is static. 

Analysis of the slope of Ti{H) at low fields (dTi/dH), 
and the exponent /? from Eq. ^as a function of p reveals 
more clearly the importance of interactions and percola- 
tion. In Fig. 21 we plot these parameters vs. p/pc- The 
dotted lines are guides to the eye. Well above Pc the 
(3 — 1, namely, the muon relaxation is a pure exponen- 
tial one, but as p decreases the nature of the muon re- 
laxation changes gradually and becomes non-exponential 
with /3 — 0.5 well below pc- Similarly dTi/dH changes 
smoothly across the percolation threshold. Neither pa- 
rameters depend critically on pc. We therefore conclude 
that while the dynamical fluctuation are sensitive to in- 
teractions, they are not governed by percolating spins 
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FIG. 3: The muon relaxation time Ti as a function of the 
longitudinal field. The solid lines are linear fits of the low 
field region at T = 100 mK. dTi/dH is obtained from the fits. 



FIG. 4; dTi/dH, and (3 vs. p/pc where pc = 0.39. The 
dotted lines are guides to the eye. 



but rather by smaller groups of spins. 

Next we present the NSE data on the same polycrys- 
talline pure {p — 1) and diluted (p — 0.37) samples as in 
the /iSR experiment. These experiments were performed 
at the Institut Laue-Langevin, France. The data from the 
pure sample only was briefly discussed in Ref. ^3 
included here, with further analysis, for completion. In 
NSE experiments, which are reviewed in Ref. jl8l |. a po- 
larized neutron beam is used. The polarization is perpen- 
dicular to the magnetic field H, which is directed along 
the neutron's momentum. The neutron spins therefore 
precess, and the difference in the number of rotations in 
two identical magnetic fields before and after the sample 
is a direct measure of a velocity change (hence, the energy 
transfer) that the neutron suffered in the scattering pro- 
cess. By varying the field one controls the Fourier time 
t = If^-fH, and the NSE method measures the scattering 
function S'(q, t) in the time domain directly, in addition 
to the spatial domain (q). The results are presented in 
Fig. panels (a) and (b) after summing over several q's 
to reduces statistical error without affecting the q reso- 
lution significantly. This is possible since the magnetic 
signal is very broad in q and NSE is a very poor q res- 
olution technique. The relaxation is found to have two 
components: A fast one, which cannot be observed by 
the NSE experiment but leads to an immediate reduc- 
tion of the normalized intermediate scattering function 
5(q, i)/S'(q, 0), and a slower one. In the pure sample 
(panel b) the relaxation of the slow component is consis- 
tent with a power law decay since it is linear on a log-log 
plot. In the diluted sample (panel a) no relaxation could 
be observed. The general evolution of the slow relaxing 
part of 5'(q, t) with doping is consistent with the fiSK 



experiment. The more diluted the sample the slower the 
fiuctuations are and the closer this fraction of the spins 
to static behavior. While this experiment does not elu- 
cidate on role of Pc it serves as a basis for a deeper un- 
derstanding of the fiSK data, especially the linear field 
dependence of the muon Ti , which is very unusual. 

In all magnets with exponentially decaying correlation 
function, Ti is a linear function of H'^ in the low field limit 
{g^sH -C J) relevant in the present case 0. Using the 
NSE results we speculate on the origin of this linear field 
dependence and its evolution with doping. We consider 
only the slow relaxing part of the NSE data, since only 
this part is relevant in the /iSR time window, and we 
make the assumption that 

^(q,i)-F(q)$(i) (2) 
where F{q) is not specified. 



X is the power of the correlation function, and Tg is an 
early time cutoff introduced to normalize the correlation 
at i = 0. However, Tg is so small that it can not be 
observed experimentally. It is therefore ignored in the 
denominator. Under this assumption the slope of the 
lines in Fig. [S] is x. The linear fit in panel (b) gives the x 
values shown in the figure. All these values are consistent 
with X — 0.02 with no q dependence. In the pure system, 
where the muon spin relaxation is exponential, Tj~^ is 
determined by the Fourier transform of the auto correla- 
tion function (S,(<)Si(0)) at the frequency / = 1^. 
This correlation function is related to the sum over q of 
S{(l,t) Therefore, under the assumption of Eq.[21we 
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FIG. 5: Neutron Spin Echo measurements in two 
(TbpYi_p)2Ti207 samples presented on a log-log scale. The 
solid lines are linear fits {x is the slope). 



have 



^{t) cos (7^ i7r) dr. 



(4) 



where 2A^ = j-^ is the RMS of the instantaneous 

transverse (to H) field distribution at the muon site (from 
the slowly fluctuating part only). This gives 



I 



2A2(2^/t,)- 



27rf 



r(l-a:)sin(— j 



In our experiment x 



Ti{H) 



0, therefore [21 
2/ 



cA2" 



(5) 



(6) 



Thus, on the basis of Eqs. 13 and 13 we can explain the 
linear field dependence of the muon Ti . We associate the 
fact that Ti > even when H = with the presence of 
static fields from other sources such nuclear moments for 
example. 



To increase confidence in Eq. |3 we check if reason- 
able values for the theoretical parameters produce the 
observed Ti. For example if Tg ^ 10"^'^ sec, Bj^ and H 
on the order of 100 G (27r/ ~ A ~ 10 MHz) and x = 0.02, 
we have (27r/Te)^ ~ 1 and Ti ~ 1 /isec, on the order of 
the measured value. We thus believe that Eq.Eldescribes 
our /iSR results appropriately. 

This allows us to conclude further that steeper Ti{H) 
implies smaller a; A^ . We argue that it is mostly x that is 
changing with decreasing p and not A. First, as observed 
in the NSE experiment, x decreases with decreasing p. 
Second, the drastic change in dTi/dH near pc cannot be 
due to a dramatic change in A. This parameter, which 
measures the local field at the muon site, certainly de- 
pends on p, but it is not sensitive to whether there is, or 
is not, percolation in the system. Therefore, the slowing 
down of spin fluctuations with decreasing p is manifested 
by a decreasing x, namely, the power law decay becomes 
closer to a constant. 

To summarize, our measurements clearly classify the 
pyrochlore (TbpYi_p)2Ti207 as a cooperative paramag- 
net. The name "paramagnet" is justified by the plateau 
in 1/Ti at low temperature. The adjective "cooperative" 
is appropriate since the spin dynamics depend on the cou- 
pling between magnetic moments. However, whether the 
spin system percolate or not is not relevant. For compar- 
ison in the kagome based compound SrCrgpGai2-9pOig 
spin dynamic was shown to be completely impartial to 
Pc E3l • We also show that the correlation of spins in this 
magnet has a power law decay with a very small power x. 
This power gets even smaller when the magnetic coverage 
decreases resulting in slowing down the spin fluctuations. 
Since the Fourier transform of behaves asymptoti- 
cally like f^~^ the fluctuation spectrum practically has 
1// behavior in the frequency range 0.5 MHz to 50 GHz. 
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